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Abstract 

Using the method of asymptotic splittings, the possible singularity structures and 

fSJ ! the corresponding asymptotic behavior of a 3-brane in a five-dimensional bulk are 

CN ■ 

^fs • classified, in the case where the bulk field content is parametrized by an analog 

lO ' of perfect fluid with an arbitrary equation of state P = ^p between the 'pressure' 

o ; 

^^ , P and the 'density' p. In this analogy with homogeneous cosmologies, the time is 

replaced by the extra coordinate transverse to the 3-brane, whose world-volume 
can have an arbitrary constant curvature. The results depend crucially on the 
constant parameter 7: (i) For 7 > —1/2, the flat brane solution suffers from a 
collapse singularity at finite distance, that disappears in the curved case, (ii) For 
7 < — 1, the singularity cannot be avoided and it becomes of the type big rip for a 
flat brane. (iii) For —1 < 7 < —1/2, the surprising result is found that while the 
curved brane solution is singular, the flat brane is not, opening the possibility for 
a revival of the self-tuning proposal. 
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1 Introduction 

In a previous work [H [2], we started a general analysis of the singularity structure and 
the corresponding asymptotic behavior of a 3-brane embedded in a five-dimensional bulk, 
using the powerful method of asymptotic splittings [3]. The main physical motivation 
is the so-called fine-tuning mechanism of the cosmological constant on the brane world- 
volume [H [S]- In Refs. in [2], we considered an extended version of the simplest model 
studied in the past, containing a massless bulk scalar field with a general coupling to 
the brane (corresponding to an arbitrary localized potential), and we allowed the brane 
world-volume to have non-vanishing constant curvature. 

We then found that the collapse singularity at a finite distance from the brane is 
present in all solutions with a fiat brane [SJ [7], but can be avoided {i.e. moved at infinite 
distance) when the brane becomes curved, either positively or negatively. The singularity 
in the fiat case is of the big bang type, characterized by the vanishing of the warp factor 
and the divergence of its derivative, as well as of the 'density' of the scalar field. 

In this paper we extend the previous analysis to the case where the bulk matter is 
described by an analog of a perfect fluid. In fact, the case of a massless bulk scalar is a 
particular case of such a fluid, corresponding to an equation of state P = 7p with 7 = 1. 
Indeed, the asymptotic behavior near the singularity of the fiat brane solution given in 
eq. (3.4) of Ref. [2], is identical to eq. (I2.18P below for 7=1. Here, we are interested 
in the general dynamics of 'evolution' of such a brane-world analog to cosmology for 
arbitrary 7, in order to reveal the various types of singularity that may develop within a 
finite distance from the original position of the brane, and on the other hand to determine 
conditions that may lead to the avoidance of the singularities shifting them at an infinite 
distance away from the brane. 

In particular, we shall show that the existence of a perfect fluid in the bulk enhances 
the dynamical possibilities of brane evolution in the fluid bulk. Such possibilities stem 
from the different possible behaviors of the fluid density and the derivative of the warp 
factor with respect to the extra dimension. The result depends crucially on the values of 



the parameter 7. In general, we find three regions of 7 leading to qualitatively different 
results: 

• In the region 7 > —1/2, the situation is very similar to the case of a massless bulk 
scalar field. A fiat brane solution has necessarily a collapse singularity at finite 
distance, which is moved at infinity when the brane becomes curved. 

• The situation changes drastically in the region — 1 < 7 < —1/2. The curved brane 
solution becomes singular while the fiat brane is regular. Thus, this region seems 
to avoid the main obstruction of the self-tuning proposal: any value of the brane 
tension is absorbed in the solution and the brane remains fiat. The main question 
is then whether there is a field theory realization of such a fluid producing naturally 
an effective equation of state of this type. 

• Finally, in the region 7 < — 1, corresponding to the analog of a phantom equation 
of state, we show that it is possible for the brane to be ripped apart in as much 
the same way as in a big rip singularity. This happens only in the flat case, while 
curved brane solutions develop a 'standard' collapse singularity. No regular solution 
is found in this region. 

Besides the above regions, the values 7 = —1/2, —1 are special: for 7 = —1/2, we find 
again a regular flat brane solution with the so-called sudden behavior |8], as well as a 
non-singular curved brane, while for 7 = —1 there is only singular curved solution. 

As mentioned above, it would be very interesting to understand whether there are 
field theory representations reproducing the 'exotic' regions of 7 < —1/2. Obviously the 
analogy of the perfect fluid concerning the positivity energy conditions does not seem to 
apply in this case where time is replaced by an additional space coordinate. However, 
some restrictions may be applied from usual field theory axioms. Also, the formation of 
singularities discussed here is better understood in a dynamical rather than the usual 
geometric sense met in general relativity. In the latter case, cosmo logical singularities 
are forming together with conjugate (or focal) points in spacetime, and for this it is 



necessary that there exists at least one timehke dimension (and any number of spacehke 
ones, greater than two). The timehke dimension forces then the geodesies to focus along 
it rather than along any of the spacelike dimensions. In the problem discussed in this 
paper, the timelike dimension is on the brane while the singularities are forming along 
spacelike dimensions in the bulk. As we show below these singularities are real in the 
dynamical sense that some component of the solution vector (a, a', p) diverges there. 
Therefore we abandon the usual interpretation according to which the universe comes to 
an end in a finite time possibly through geodesic refocusing, and instead we study how 
dynamical effects guide our brane systems to extreme behaviors. 

The structure of this paper is the following. In Section 2, we first choose appropriate 
variables and rewrite the basic field equations of the problem in the form of a dynamical 
system; secondly, we introduce some convenient terminology for the different types of 
singularity to be met later in our analysis; thirdly, we single out the possible dominant 
balances, organizing centers of all the evolutionary behaviors that fully characterize 
our problem. In Sections 3 and 4, we study carefully the asymptotics around collapse 
singularities of two types, that we call I and II, respectively. In Section 5, we explore 
the dynamics as the brane approaches a big rip singularity, while in Section 6 we look 
at a milder singularity that resembles in many ways the so-called sudden (non-singular) 
behavior introduced in Ref. [Sj. In Section 7, we analyze the possibility of avoiding 
finite- distance singularities leading to the existence of regular brane evolution in the 
bulk and finally, in Section 8 we conclude and refer to future work. 

2 Dynamics in a perfect fluid bulk 

In this Section we rewrite the brane model living in a bulk filled with a perfect fluid as a 
dynamical system in three basic variables and completely identify the principal modes of 
approach to its singularities, that is we find all the dominant balances of the system. We 
consider a three-brane embedded in a five- dimensional bulk space filled with a perfect 



fluid with equation of state P = 7p, wliere tlie pressure P and tlie density p are functions 
only of the fifth dimension, denoted by Y. We assume a bulk metric of the form 

g, = a\Y)g, + dY\ (2.1) 

where g^ is the four-dimensional flat, de Sitter or anti de Sitter metric, i.e., 

94 = -de + flg^, (2.2) 

where 

(73 = dr^ + hlg^ (2.3) 

and 

g^ = de^ + sin^ ed^"^. (2.4) 

Here /« = r, sinr, sinhr, and h^ = 1 , cosh.{Ht) / H , cos{Ht) / H {H^^ is the de Sitter 
curvature radius). We also assume an energy- momentum tensor of the form Tj^b = 
{p + P)uaUb — PgAB, where A,B = 1,2,3,4,5 and Ua = (0,0,0,0,1), with the 5th 
coordinate corresponding to Y. 

The five-dimensional Einstein equations. 

Gab = 4Tab, (2.5) 

with Kg = Mg"^ and M5 the five dimensional Planck mass, can be written in the following 
form: 

a o 

— = -i^P+^, (2-7) 

where k = 0,±1, and the prime (') denotes differentiation with respect to Y. The 
equation of conservation, 

VbT^^ = 0, (2.8) 

becomes, 

p' + 4(1 + 7)iJp = 0. (2.9) 



Introducing the new variables 



(2.10) 



Eqs. ([21]) and ([23]) take the form 



X = y, (2.11) 

.' = -2Aii±M,,, (2.12) 

/ = -4(1+7)^2, (2.13) 

X 



while eq. (12. 7p reads 



^ = -A. + --, A = 4/A. (2.14) 



y' 2 , kH' 

x^ 3 x^ 

Since this last equation does not contain derivatives with respect to Y , it is a velocity 
independent constraint equation for the system (I2.1ip - (l2.13p . Before we proceed with 
the analysis of the above system, we introduce the following terminology for the possible 
singularities to occur at a finite-distance from the brane. Specifically, we call a state 
where: 

i) a — )■ 0, a' — )■ oo and p — )■ oo: a singularity of collapse type I. 

ii) a — )■ 0, a' — )■ a'^ and p — ;■ 0: a singularity of collapse type Ila, 
a — 7- 0, a' — )■ a'g and p ^ ps'. a singularity of collapse type lib, 
a — > 0, a' — )■ a'^ and p — > oo: a singularity of collapse type lie, 
where a'^, ps are non- vanishing finite constants. 

iii) a — 7- oo, a' — 7- — oo and p — )■ oo: a big rip singularity. 

We note that with this terminology the finite-distance singularity studied in [2] is a 
singularity of collapse type I. 

The next step is to apply the method of asymptotic splittings in an effort to find all 
possible asymptotic behaviors of the dynamical system (I2.1ip -( l2.13p with the constraint 



fl2.14p . by building series expansions of the solutions around the presumed position of 
the singularity at Ys. 

We note that the system fl2.1ip - (l2.13p is a weight homogeneous system determined 
by the vector field 

f = (^^,-2A^i±Mzx,-4(l + 7)^^y. (2.15) 

In order to compute all possible dominant balances that describe the principal asymp- 
totics of the system we look for pairs of the form, 

i3 = {a, p}, where a={a,f3,S), P={p,(l,r), (2.16) 

with 

{p,q,r)EQ^ and {a, (3,5) E C \ {0}, (2.17) 

by setting {x,y,z) = {aTP, PT'^,6T') in the system (l2llD - (l27[3l) . where T = F - F, is 
the distance from the singularity. We find after some calculation the following list of all 
possible balances for our basic system fl2.1ip -( l2.14p : 

7^1 = |(a,ap,^p''),(p,p-l, -2)}, ^=^^^^^,77^ -1/2,-1,(2.18) 

^B2 = {(a,a,0),(l,0,-2)}, 7 7^-1/2, (2.19) 

-i/2i33 = {(a,a,0),(l,0,r)}, (2.20) 

-1/2^4 = {(«,«, (5), (1,0, -2)}, (2.21) 

-1/2^5 = {(«,0,0),(0,-l,r)}, (2.22) 

where _i/2'Bj =^=_i/2 ^i- Notice that, as already mentioned in the introduction, the first 
balance ^Bi for 7 = 1 coincides with the one found in [2] in eq. (3.4), where the fluid 
was replaced by a massless bulk scalar field with an arbitrary coupling to the brane. 

The above balances are exact solutions of the system and they must therefore also 
satisfy the constraint equation f l2.14p . This fact alters the presumed generality of the 
solution represented by each one of the balances above and determines uniquely the 
type of spatial geometry that we must consider: The balances ^Bi and -1/2-65 are found 



when we set fc = 0, and describe a (potentially general) solution corresponding to a flat 
brane, while the balances ^82 and _i/2'B3 were found when k ^ and describe particular 
solutions of curved branes (since we already have to sacrifice the arbitrary constant a by 
imposing a^ = kH"^). For the balance _i/2'B4, on the other hand, k is not specified and 
hence it describes a particular solution for a curved or fiat brane (particularly since we 
have to set 5 = (3/(2v4))(l - kH^/a"^) to satisfy eq. (EH])). 

Each one of these balances are analyzed in detail in the following sections according 
to the nature of asymptotic behaviors they imply. 

3 Collapse type I singularity 

We shall focus in this Section exclusively on the balance ^Bi and show that for certain 
ranges of 7 it gives the generic asymptotic behavior of a flat brane to a singularity of 
collapse type I. Our analysis implies that such behavior can only result from a -yBi type 
of balance. 

Our purpose is to construct asymptotic expansions of solutions to the dynamical 
system (I2.1ip -( l2.14p in the form of a series solution defined by 

x = TP(a + E~iC,T^>), (3.1) 

where x = (x, y, z), Cj = (cji, Cj2, Cjs), and s is in this case the least common multiple of 
the denominators of the positive /C-exponents (cf. [3], [9]). As a first step we calculate 
the Kowalevskaya matrix, /C = Df (a) — diag(p), where -Df (a) is the Jacobian matrix of 



•B 



fH: We have, 



/ 



Df{x,y,z) 



\ 



--{l + 2^)Az 



--(l + 27)Ax 



(3.2) 



V ^^'+^^1 -'^^+^)^ -'(^+^)! y 

to be evaluated on a. The balance ^Bi has a = (a, ap, 3p'^/2A), and p = {p,p — 1, —2), 
with p = 1/(2(7 + !))• Thus the Kowalevskaya matrix (/C-niatrix in short) for this 
balance, is 



y/Ci = Di [a, ap, — / ) - diag(p, p-1,-2) 



Di a 



2(l + 7)'8A(l + 7)^ 



— diag 



2(1 + 7) 



1 + 27 



1 1 + 27 

2(1 + 7)' 2(1 + 7)' 

^ 



-2 



1 + 27 

'4(1+7)2 2(1+7) 



:i + 2-f)Aa 



\ 4(l + 7)Ma 2(l + 7)Aa 



(3.3) 



Let us now calculate what the /C-exponents for this balance actually are. Recall 
that these exponents are the eigenvalues of the matrix ^/Ci and constitute its spectrum, 
spec{.y}Ci). The arbitrary constants of any (particular or general) solution first appear 
in those terms whose coefficients c^ have indices k = gs, where g is a non-negative 
/C-exponent. The number of non-negative /C-exponents equals therefore the number of 
arbitrary constants that appear in the series expansions of (13. ip . There is always the —1 



^f is the vector field resulting froni the dynamical system (|2.1ip - (|2.14l) and {a, p} is the balance ^Bi. 



exponent that corresponds to an arbitrary constant that is the position of the singularity 
at Yg. The balance ^i3i corresponds thus to a general solution in our case if and only if 
it possesses two non-negative /C-exponents (the third arbitrary constant is the position 
of the singularity, Ys). Here we find 

spec(,/Ci)=|-l,0,^^^|. (3.4) 

The last eigenvalue is a function of the 7 parameter and it is positive when either 7 < — 1 , 
or 7 > —1/2. We consider here the case 7 > —1/2 since, as it will soon follow, this range 
of 7 is adequate for the occurrence of a collapse type I singularity. The case of 7 < — 1 
leads to a big rip singularity and will be examined in Section 5. 
Let us assume 7 = —1/4 for concreteness. Then 

.1/4^1 = {(«,-2«/3,2/(3A)),(2/3,-l/3,-2)}, (3.5) 

spec(_i/4/Ci) = {-1,0,2/3}. (3.6) 

Substituting in the system (12. lip - ( 12713]) the particular value 7 = —1/4 and the forms 

X = E°lo9lT^■/'+'/^ y = S°^oC,2T^■/^-^/^ z = S°loC,3T^■/^-^ (3.7) 

we arrive at the following asymptotic expansions: 

X = aT'/' - ^c,,T'/' + ■ ■ ■ , (3.8) 

y = ^«T-i/3_2^«c23Ti/=^ + ---, (3.9) 

z = ^T-' + C23T'^^' + --- . (3.10) 

For this to be a valid solution we need to check whether for each j satisfying j/3 = g 
with g a positive eigenvalue, the corresponding eigenvector v of the _i/4/Ci matrix is such 
that the compatibility conditions hold, namely, we must have 



v'^-Pj = 0, (3.11) 



where Pj are polynomials in Cj, . . . , Cj_i given by 

-1/4IC1CJ - {j/3)cj = Pj. 



(3.12) 



Here the corresponding relation j/3 = 2/3, is valid only for j = 2 and the compatibility 
condition indeed holds since, 



_l/4/Ci - (2/3)X3)C2 



/ 4 



2 

'9 



1 
3 



\ 3Aa Aa 
Eqs. dSSD-dSin]) then imply that as T ^ 0, 



^ 

Aa 

2 
"3 / 



C23 



2 ' ^ 

2Aa 



(3.13) 



voy 



a — )■ 0, a' — )■ 00, p — )■ 00. 



(3.14) 



This asymptotic behavior corresponds to a general solution of a flat brane that is valid 
around a collapse I singularity. We thus regain a behavior similar to the one met in [2] 
for the case of a flat brane in a scalar field bulk. 

4 Collapse type II singularities 

In this Section, we show that for a curved brane (A; = ±1) the long-term (distance) be- 
havior of all solutions which depend on the asymptotics near finite-distance singularities 
turn out to be of a very different nature. In particular, we shall show that the balances 
^82 for 7 < —1/2, _i/2i33 for r < — 2 and _i/2'B4 (as we have already mentioned ^82 
and -1/2B3 correspond to a curved brane whereas -1/2-64 corresponds to a flat or curved 
brane), imply the existence of a collapse type Ila, b or c singularity. This is in sharp 
contrast to the asymptotic behavior found for a curved brane in the presence of a bulk 
scalar field (cf. ||2jj), wherein there are no finite-distance singularities. 



10 



For the balance ■yB2 we find that 



y/Ca = Df (a, a, 0) - diag (1, 0, -2) 



/ 



\ 



-1 1 

-^Aa(l + 27) 

\^ -2(1 + 27) J 



(4.1) 



and hence, 

spec(^/C2) = {-1, 0, -2(1 + 27)}. (4.2) 

We note that the third arbitrary constant appears at the value j = —2(1 + 27), 7 < ~ V^- 
After substituting the forms. 



X 






(4.3) 



in the system (I2.1ip - (l2.13p . to proceed we may try giving different values to 7: Inserting 
the value 7 = —3/4 in the system for concreteness we meet a third arbitrary constant at 
j = I (spec(_3/4/C2) = { — 1,0, 1}). We then arrive at the following asymptotic forms of 
the solution: 



X = ai + -—C13T + 

D 

y = a + — CigTH 

z = CisT"^ H , 



(4.4) 

(4.5) 
(4.6) 



iB 



where Ci 3 7^ q We need to check the validity of the compatibility condition for j 
But this is trivially satisfied since, 

/ -2 1 \ 



-3/4 



/C2 - X3)Ci 



-1 Aa/3 
\ / 



/ Aa/6 \ 



Cl3 



V 



Aa/3 
1 



/o\ 





(4.7) 



^If we do not set from the beginning 7 = —3/4 but instead we let 7 be arbitrary, then in the last 
step of the calculations at the j — I level we find that either C13 = or 7 = —3/4. 
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The series expansions in eqs. fl4.4l) - fl4.6p are therefore vahd and we conclude that as 

a — !■ 0, a' — )■ a, p — )• oo, a 7^ 0. (4.8) 

This is a collapse type lie singularity. It will follow from the analysis below that the 
behavior of p depends on our choice of 7 (thus giving rise to three possible subcases of 
a type II singularity). Indeed, choosing for instance 7 = — 1 (spec(_i/C2) = { — 1!0,2}), 
we find that the solution is given by the forms, 



X 

y 

z 



al+ -^C23T + 
y 

« + ^C23T +■■ 
C23H , 



(4.9) 

(4.10) 
(4.11) 



,B 



where C23 7^ 0|j. Note that the compatibility condition is satisfied here as well since. 



/ 



:_i/c2 - 2x3)02 



We see that as T — )■ 0, 



\ 



-2 2Aa/3 
\ / 



C23 



/ Aa/9 \ 

Aa/3 

V 1 / 



/o\ 





(4.12) 



a — )■ 0, a' — !■ a, p — )■ C23, a 7^ 0. 



(4.13) 



This is a collapse type lib singularity in our terminology and is clearly different from 

(iSD. 

A yet different behavior is met if we choose for instance 7 = —5/4. The /C-exponents 
are given by spec(_5/4/C2) = {—1)0,3}, and the series expansions become, 

Aa 



X 

y 

z 



Aa ^„o 
a + — C33T=^ + -- 

c^{^ H , . 



(4.14) 

(4.15) 
(4.16) 



■^Had we let 7 be arbitrary we would have found that in the step j = 2 of the procedure either C2 3 = 
or 7 = — 1. 
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where C33 7^ Q These expansions are vahd locally around the singularity since the 
compatibility condition holds true because, 

/ -4 1 \ / AajVl \ / \ 

(-5/4/C2 - 3X3)01 = -3 Aa C32 AajZ 

V ; V 1 / 

For T — )■ 0, we have that 







(4.17) 



a — )■ 0, a' — )■ a, p — ;■ 0, a 7^ 0, 



(4.18) 



which means that this is a collapse type Ila singularity. This balance therefore leads to 
the asymptotic behavior of a particular solution describing a curved brane approaching 
a collapse type II singularity, i.e., a — ?■ and a' — )■ a. The behavior of the density of 
the perfect fluid varies dramatically: we can have an infinite density, a constant density, 
or even no flow of "energy" at all as we approach the flnite-distance singularity into the 
extra dimension at Ys-, depending on the values of the 7 parameter. 

We now turn to an analysis of the balances _i/2i33, for r < —2, and _i/2'B4- The 
/C-matrix for -1/2-63 is 



/ 



-i/2A^3 = -Df («, «, 0) - diag(l, 0, r) 



4 1 






\ 







2- 


t) 



(4.19) 




\ 

and hence, 

spec(_i/2/C3) = {-1, 0, -2 - r}. 

Taking —2 — r > 0, we have two non-negative /C-exponents. (The case —2 — r < is 
considered later, in Section 8, since it is quite different, it does not imply the existence 
of a flnite-distance singularity.) For r = —3 as an example, we substitute the forms 



(4.20) 



X 






y 



Sr=oC,2T^ 



z 



sr=o93T^-^ 



(4.21) 



^Here again if had let 7 be arbitrary we would have found that in the step j = 3 of the procedure 
either C33 = or 7 = —5/4. 
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and arrive at the expansions 



X 

y 

z 



aT + ■ ■ • . 
a^ , 



(4.22) 
(4.23) 
(4.24) 



The compatibihty condition is satisfied because 



/ 



1/2^3 - '^z)Cl 



1 o\ 



/o\ 



Cl3 







/o\ 







(4.25) 



0-10 
\ 0/ 

and so the expansions (I4.22p - (l4.24p are vahd ones in the vicinity of the singularity. The 
general behavior of the solution is then characterized by the asymptotic forms 



a — )• 0, a' — > a, p — )• oo, a 7^ 0. 



(4.26) 



The balance _i/2'B3 for r < —2 implies therefore the existence of a collapse type lie 
singularity during the dynamical evolution of the curved brane living (and moving) in 
this specific perfect fiuid bulk. 

The balance -i/2'B4 on the other hand is one with 

/-I 

-1/2/C4 = -Df (a, a, 5) — diag(l, 0, —2) = 







\ 



v 





a 





-H^ 

a 



(4.27) 



/ 



and 



spec(_i/2/C4) = {-1,0,0}. 



(4.28) 



We note that the double multiplicity of the zero eigenvalue refiects the fact that there 
were already two arbitrary constants, a and 5 in this balance (recall though that 5 had 
to be sacrificed in order for this balance to satisfy the constraint (I2.14p ). We can thus 
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X = 


= aT + --- 


y = 


= « + ■■■ , 


Z — 


= 5T~^ + - 



write 

(4.29) 
(4.30) 
(4.31) 

so that as T — )■ 0, a collapse type lie singularity develops, i.e., 

a -)■ 0, a' -)■ a, p -)■ oo, a 7^ 0. (4.32) 

5 Big rip singularities 

In this Section we return to the balance ^Bi but focus on different 7 values. In particular, 
we show that when 7 < —1, a flat brane develops a big rip singularity in a finite distance. 
This new asymptotic behavior implied by the balance ^Bi (when 7 < —1) is equally 
general to the one found in Section 3. 

For purposes of illustration, let us take 7 = —2. Then the balance -2B1 and the 
_ 2 /Ci -exponents read, respectively, 

_2-Bi = {(«,-«/2,3/(8A)),(-l/2,-3/2,-2)}, (5.1) 

spec(_2/Ci) = {-1,0,3}. (5.2) 

Substituting the value 7 = — 2 in our basic system given by eqs. (I2.11l) - (l2.13p . and also 
the forms 



X 



^T=o^ji^'~'^\ y = E- oC,2T^-'/^ ^ = E~ o93T^■-^ (5.3) 



we expect to meet the third arbitrary constant at j = 3. Indeed we find: 

X = aT''/^ + lAacssr'/^ + --- , (5.4) 

y = -|T-^/^ + ^Aac33T3/^ + ---, (5.5) 

Z = ^T-2^C33T + ---, C33^0. (5.6) 
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The compatibility condition is trivially satisfied for j = 3, since the product (-2^1 
3X3)03 is identically zero: 



/ -^ 



(_2/Ci - 3X3)03 



\ n-Aa\ 
3 



3 3 

- — 2Aa 

4 2 



v 



C33 



/ 



-Aa 
3 



V 1 / 







(5.7) 



iAa 2Aa 

The series expansions given by eqs. fl5.4p - (15.61) are therefore valid asymptotically for 
T — )■ so that we end up with the asymptotic forms 



a — )■ 00, a — )■ —00, p — )■ cxD. 



(5.^ 



We therefore conclude that the balance ^Bi leads to a general solution in which a flat 
brane develops a big rip singularity after 'traveling' for a finite distance when the bulk 
perfect fluid satisfies a phantom-like equation of state, i.e., 7 < — 1. Note that using the 
analogy between the warp factor of our braneworld and the scale factor of an expanding 
universe, we can say that this singularity bares many similarities to the one studied 
in Refs. [10], [11], [12], since it is also characterized by all quantities a, a', p, and 
consequently P, becoming asymptotically divergent. Thus, the results in this Section 
indicate that a fiat brane traveling in a 7 < —1 fluid bulk develops a big rip singularity. 
This implements the behavior found in Section 3 of the present paper, wherein the same 
brane moving in a 7 > —1/2 fluid bulk 'disappears' in a big bang-type singularity. 
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6 Sudden behavior 

As our penultimate mode of approach to the finite-distance singularity, we examine here 
the balance -1/2B5 = {{a, 0, 0), (0, —1, r)}. This balance has 

/ 1 \ 



.i/2/C5 = /^f(a,0,0)-diag(0,-l,r) 



1 
\ -r / 



(6.1) 



and 



(6.2) 



spec(_i/2/C5) = {1,0, -r}, 

so we shall have to set r = 1 in order to have the necessary —1 eigenvalue corresponding 
to the arbitrary position of the "singularity", Yg. After substitution of the forms 






we find that the solution reads 



(6.3) 



X = a + CiiT + 

y = ci 1 H , 

z = + --- . 



(6.4) 
(6.5) 
(6.6) 



The compatibility condition is satisfied since 



^1/2^5 ~ 2^3)01 



and we see that as T — )• 0, 



/-I 1 o\ 




/l\ 




/o\ 





Cll 


1 


= 





\ yi 




UJ 




U/ 



X — )• ct, y ^ Cll, -z — )■ 0, 0:7^0. 



(6.7) 



(6. 



This clearly indicates that the brane experiences the so-called sudden behavior (cf. [8]) 
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7 Behavior at infinity 

A totally different picture than what we have already encountered up to now in our 
analysis of brane singularities in a fluid bulk, is attained using the balance ^Bi with 
— 1 < 7 < —1/2, or the balance ^82 with 7 > —1/2, or the balance _i/2'B3 with r > —2. 
We show in this section that these three balances and only these offer the possibility 
of avoiding the finite-distance singularities met before and may describe the behavior of 
our model at infinity. 

We begin with the balance jBi when —1 < 7 < —1/2. Choosing for instance 
7 = —4/5, we find spec(_4/5/Ci) = {-1,0,-3} and hence we may expand {x,y,z) in 
descending powers in order to meet the arbitrary constant appearing at j = —1 and 
j = -3, i.e.. 



X 






y 









YJ~2 



(7.1) 



We find: 



X = aT'/' + c^^^T'/' + 3/{10a)c\,T'/' + c^slr-'/' + ■■■, (7.2) 

y = 5a/2T3/2^3/2c_iiT^/2 + 3/(20a)c'_iiT-^/'-l/2c_3iT-'/=' + --- , (7.3) 
z = 75/{8A)T-^-15/{2Aa)c-iiT-^ + 9/{2Aa^)c\^T'^ + 

+ (-15/(2Aa)c_3i - 9/iAAa^)ci,,) T'^ + ■ ■ ■ . (7.4) 



The compatibility conditions at j = — 1 is satisfied since 



-4/5^1+2^3)C-l 



/ -3/2 


1 \ 




( ^ ] 




(o\ 


15/4 


-1/2 2Aa/5 


C-ll 


3/2 


= 





\ 75/(4Aa) 


-15/{2Aa) 1 ) 




i^ -15/{2Aa) ) 




(7.5) 
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But for j = — 3 we find 



-4/5^^1 + 3X3)C_3 = 

/ 1/2 1 ^ \ ^ 

15/4 3/2 2Aa/5 

\ 75/(4Aa) -15/(2Aa) 3 y 

( ° ^ 

-9/(10a2)c3_ii =P-3. 

\ -27/ (A Aa^)c^_^^ J 



C-31 \ 

-l/2c_3i 



(7.6) 



An eigenvector corresponding to the eigenvalue j = —3 is v^ = (— 2Aa/15, Aa/lS, 1), 
and hence we have 

v^ ■ P-3 ^ 0, (7.7) 

unless c_ii = 0. In order to satisfy the compatibility condition at j = —3 we set 
c_ii = 0. The solution ([73)-(I73D with c_ii = reads 



y = 5a/2T3/2_^/2c_3iT-3/2 + -- 



z 



75/(8A)T-2 - 15/(2Aa)c_3iT-^ + 



(7.8) 

(7.9) 

(7.10) 



and it is a particular solution containing two arbitrary constants. As S = 1/T — ;■ oo, we 
conclude that 



a — >■ oo, a — )■ oo, p — )• oo. 



(7.11) 



and we can therefore avoid the finite-distance singularity in this case. 

Next we examine the balance ■./B2 when 7 > —1/2. For 7 = 0, we have that 
spec(o/C2) = { — 1,0, —2}, and hence we substitute 



X 






(7.12) 
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and find: 



X = aT + c^i 1 — Aa/3c_23T' "^ + ■ • • , 
y = a + Aa/ 3c -23'^"'^ -\ , 

The compatibility conditions at j = —1 and j = —2 are indeed satisfied since 

/O 1 \ /i\ /o\ 



(7.13) 
(7.14) 
(7.15) 



]/C2 +X3)c_i 



1 -2Aa/3 

Vo -1 / 



c-11 







and 



/l 1 \ 



(0/C2 + 2X3)C_2 



/ -Aa/3 \ 



2 -2Aa/3 
\ / 

As S* = 1/T — )> 00, we conclude that 



C-2 3 



V 



Aa/3 
1 





VO/ 


VO/ 



(7.16) 



(7.17) 



a — > 00, a — )■ 00, p — 7- 00, 



(7.18) 



and the finite-distance singularity in shifted at an infinite distance. 

We now move on to the balance _i/2'B3, r > —2. In this case this balance has two 
negative /C-exponents. If we choose the value r = 0, then the spectrum is found to be 



spec(_i/2/C3) = {-1,0,-2}, 



and so inserting the forms 



X 






we obtain 



X = aT + c_ii, 
y = a, 

z 



C-23T-2 + 



(7.19) 

(7.20) 

(7.21) 
(7.22) 
(7.23) 
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which vahdates the coinpatibihty conditions at j = — 1 and j 

/O 1 \ /l\ 



1^-1/2^3 +2r3)c_i 



and 



(_l/2/C3 + 2I3)C_2 

We see that as S = 1/T — ;■ oo, 



1 

V -1 / 

/l 1 o\ 

2 
\ y 



C-ll 



C-2 3 





/o\ 





-2 since 

/o\ 

, 



/o\ 



vo/ 



(7.24) 



(7.25) 



a— )'C_ii, a' — )■ a, p — )■ oo, 07^0, 



(7.26) 



so that the balance -1/2-63 for r > —2 also offers the possibility of escaping the finite- 
distance singularities. Hence in such cases we find a regular (singularity free) evolution 
of the brane as it travels in the bulk filled with the type of matter considered above. 

8 Conclusions 

We have studied the dynamical 'evolution' of a braneworld that consists of a three-brane 
embedded in a five- dimensional bulk spacetime filled with a 'perfect fluid' possessing a 
general equation of state P = 7p, characterized by the constant parameter 7. 

For a flat brane we find that it is possible to have within finite distance from the 
brane the analogous type of collapse singularity met previously in [2J. We call this a 
collapse type I singularity and describe its nature using the behavior of the warp factor, 
its derivative and the density of the fluid, asa— )-0, a'— 7'00,p— >oo. In|2], the bulk 
matter was modeled by a scalar field and this singularity was the only type possible. 
Here we showed that when the scalar field is replaced with a perfect fluid, in addition 
to that singularity which appears inevitably in all flat brane solutions with 7 > —1/2, 
there are two other new types (for a flat brane): The first one is the very distinct big 
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rip singularity which occurs with a — )■ oo, a' — )■ — oo, p — )■ oo and only when a phantom 
type equation of state with 7 < — 1 is considered. The second one is a collapse type 
lie singularity which may be described by the behavior a — )■ 0, a' — )■ a and p — )■ 00. 
We note that this latter singularity is less general than the collapse type I and the big 
rip singularities and it arises only when 7 = —1/2. Besides these singular solutions, 
we found the surprising result of flat branes without finite-distance singularities in the 
region —1 < 7 < —1/2. Moreover, for 7 = —1/2 the solution has the sudden behavior 
with a and a' finite and vanishing density p — )■ [H]. 

In contrast to the bulk scalar field case where all curved brane solutions were regular, 
now we found also singular such solutions. The possible corresponding finite-distance 
singularities are the ones comprising the collapse type II class. These are singularities 
with a — > 0, a' — )■ a and p —t- 0, ps, 00 (corresponding to types Ila, b and c respectively). 
The interesting feature of this class of singularities is that it allows the 'energy' leak into 
the extra dimension to vary and be monitored each time by the 7 parameter that defines 
the type of fluid; they all arise in the region 7 < —1/2. On the other hand, we showed 
that for a curved brane the possibility of avoiding the finite-distance singularities that 
was offered in [2] is still valid here, but only in the region 7 > —1/2. 

For illustration, we present a summary of all different behaviors we found for flat 
and curved branes in the table below, using the notation for the various singularities 
introduced in Section 2 after eq. fl2.14p and the balances (I2.18p -f l2^22|) . 
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equation of state 


flat brane 


curved brane 


P = 1P 


type 


balance 


type 


balance 


7>-l/2 


singular type I 


,Bi 


regular 


^^2 at oo 


7 = -1/2 


singular lie 
regular sudden 


-1/2^4 
-1/2^5 


regular 
singular lie 


-1/2-63 , r > -2 at cx) 
-1/2-63 , r < -2 or -1/2-64 


-1 <7< -1/2 


regular 


^Bi at oo 


singular IIa,b,c 


,B2 


7 = -l 


no solution 




singular lib 


,B2 


7<-l 


singular big rip 


,Bi 


singular Ila 


,B2 



An open question is whether there exist physical constraints on 7 analog to the 
weak and strong energy conditions of matter perfect fluid in ordinary cosmology. A 
related question is to find possible field theory realizations of the 'exotic' regions of 
7 < —1/2, where interesting solutions with unexpected behavior were found. The most 
important issue of course is to clarify the possibility of singularity avoidance at finite 
distance in flat brane solutions. There is no reason why the non-singular behavior for 
flat branes discovered here should not persist for arbitrary values of the brane tension 
and, indeed, it is to be expected that only particular asymptotic modes of behavior, 
that is specific detailed forms of asymptotic solutions, would depend on such values. 
Thus, the self-tuning mechanism appears to be a property of a general (non-singular) 
flat brane solution, that depends on two arbitrary constants in the region — 1<7<— 1/2 
(three for the general solution with sudden behavior when 7 = —1/2). Similarly, as we 
have shown here, the existence of a singularity may be independent of the sign of the 
scalar curvature (as long as the latter remains non zero for curved branes), but the 
particular way of asymptotic approach to the singularity is sensitive to that sign and it 
may therefore change with different values of brane tension. 

It would also be interesting to further investigate whether the properties of finite- 
distance singularities (and their possible avoidance) encountered here continue to emerge 
in more general systems, such as the case in which a scalar field coexists with a perfect 
fluid in the bulk [13] . The analysis of this more involved case may also shed light to the 
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factors that control how these two bulk matter components compete on approach to the 
singularity, or even predict new types of singularities that might then become feasible, 
as well as possible situations where they can be avoided. 
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